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SPECTRAL EVOLUTION OF THE SEPTEMBER 2010 GAMMA-RAY FLARE FROM THE CRAB NEBULA 
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ABSTRACT 

Strong gamma-ray flares from the Crab Nebula have been recently discovered by AGILE and con- 
firmed by Fcrmi-LAT. We study here the spectral evolution in the gamma-ray energy range above 
50 MeV of the September 2010 flare that was simultaneously detected by AGILE and Fermi-LAT. 
We revisit the AGILE spectral data, and present an emission model based on rapid (within 1 day) 
acceleration followed by synchrotron cooling. We show that this model successfully explains both the 
published AGILE and Fermi-LAT spectral data showing a rapid rise and a decay within 2-3 days. 
Our analysis constrains the acceleration timescale and mechanism, the local magnetic field, and the 
particle distribution function. The combination of very rapid acceleration, emission well above 100 
MeV, and the spectral evolution consistent with synchrotron cooling contradicts the idealized scenario 
predicting an exponential cutoff at photon energies above 100 MeV. We also consider a variation of 
our model based on even shorter acceleration and decay timescales which can be consistent with the 
published averaged properties. 
Subject headings: 



1. INTRODUCTION 

The Crab Nebula is at the center of the SN1054 su- 
pernova remnant and consists of a rotationally-powered 
pulsar interacting with a surrounding nebula through a 
relativistic particle wind (e.g., Hester 2008). The Crab 
pulsar is quite powerful (of spindown luminosity Lpsr= 
5-10 38 erg s _1 , and spin period P — 33 ms), and is en- 
ergizing the whole nebula with its wave/particle out- 
put. The inner nebula shows distinctive optical and 
X-ray brightness enhancements ("wisps", "knots", and 
the "anvil" aligned with the pulsar "jet") (Scargle 1969; 
Hester 1995, 2002, 2008; Weisskopf 2000). These local 
variations have been attributed to enhancements of the 
synchrotron emission produced by instabilities and/or 
shocks in the pulsar wind outflow. However, when av- 
eraged over the whole inner region (several arcminutc 
across) the Crab Nebula has been considered essen- 
tially stable, and used as a "standard candle" in high- 
energy astrophysics. The Crab Nebula X-ray contin- 
uum and gamma-rays up to ~ 100 MeV energies are 
modelled by synchrotron radiation of accelerated parti- 
cles in an average nebular magnetic field B = 200 /iG 
(Hester 2002, deJager et al. 1996, Atoyan & Aharo- 
nian 1996, Meyer et al. 2010). Emission from GeV to 
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TeV energies is interpreted as inverse Compton radiation 
by electrons/positrons scattering CMB and nebular soft 
photons (deJager & Harding 1992, deJager etal. 1996, 
Atoyan & Aharonian 1996, Meyer et al. 2010). 

Decades of theoretical modelling of this system 
(e.g., (Rees & Gunn 1974, Kennel & Coroniti 1984, 
deJager & Harding 1992, deJager etal. 1996, Atoyan 
& Aharonian 1996, Arons 2008, Meyer et al. 2010) 
offer the picture of a remarkable nebular system 
energized by a MHD pulsar wind interacting with 
the environment through a sequence of "shocks" or 
dissipation features localized at distances larger than 
a few times 10 17 cm. Efficient particle acceleration at 
the pulsar wind termination shock regions is believed 
to be occurring either thro ug h diffusive processes, e.g. , 
dBlandford fc Eichler 19871 : Ide Jager fc Harding 1993: 
Ide Jager et al. 19961 : lAtovan fc Aharonian 1996f l . 

shock-drift acceleration (e.g., Begelman & 
Kirk 1990) or ion- mediated acceleration (e.g., 
(|Arons 20081; ISpitkovskv fe Arons 2001 V Several 
diffusive acceleration models imply acceleration rates 
of orde r of the relativistic electron cyclotron frequency 
(e.g.. (Ide Jager fc Harding 1992t Ide Jager et al. 1996t 
lAtovan fc Aharonian 1996T )L Assuming equality be- 
tween the accelerating electric field and the magnetic 
field at the acceleration site and synchrotron cooling 
in the co-spatial magnetic field leads to a most cited 
constraint for the maximum radiated photon energy 
(e.g., Aharonian 2004) 
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with a = e 2 /he the fine structure constant, c the speed of 
light, and m e the electron's mass. Eq. [T]applies in a nat- 
ural way to diffusively accelerated particles, and E 1%rnax 
turns out to be independent of the local magnetic field. 
According to the assumptions underlying this formula, 
emission above 100 MeV would be difficult to sustain 
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in the Crab Nebula environment. Indeed, the exponen- 
tial cutoff shown by the average gamma-ray spectrum in 
the 10 MeV - 10 GeV range supports this ide alization 
(|de Jaeer fc Harding 19921: Ide Jager et al. 19961) . 

However, the recent discovery by the AGILE satel- 
lite of a strong gamma-ray flare ab ove 100 MeV from 
the Crab Nebula in September 2010 (|Tavani et al. 20101: 
iTavani et al. 20fl|) and th e confirmation by t he Fermi- 
LAT (jBuehler et al. 20101: lAbdo et al. 201 ID substan- 
tially change this picture. Three substantial gamma-ray 
flaring episodes from the Crab Nebula have been an- 
nounced so far (Tavani et al. 2011, hereafter Til; Abdo 
et al. 2011, hereafter All). The flaring activity was de- 
tected only in the gamma-ray energy range 100 MeV - 
a few GeV, and it is attributed to transient nebular un- 
pulsed emission. No significant variations were detected 
in other bands, except for the local enhancement in the 
"anvil" region revealed by high spatial resolution obser- 
vations of HST and chandra. 

Three features of the September 2010 event are rele- 
vant: (1) the event develops within 3-4 days (whereas 
the others last about 2 weeks); (2) the gamma-ray rise- 
time appears to be remarkably short, r < 1 day (Til); 
(3) the flaring gamma-ray spectrum extends well above 
the limit of Eq. 1 (Til, All). A flare production site in 
the inner nebula of size L < 10 16 cm is favored by both 
the peak isotropic gamma-ray luminosity L p k, 5T0 35 
erg s _1 (which implies for a (3-5) % radiation efficiency 
that about (2-3)% of the total spindown pulsar lumi- 
nosity is dissipated at the flaring site) and by the flare 
risetime of ~1 day. We noticed that the "anvil region" 
P'knot -2" and possibly "k not-1") in the Crab Nebula 
dScargle 1 969: Hester 2008) is an excellent flare site can- 
didate also because of its alignment with the relativistic 
pulsar jet (Til). 

A number of important theoretical questions are raised 
by these detections. However, the published spectra 
of the September 2010 event are not homogeneous be- 
cause of different integration times: a 2-day timescale 
for the AGILE data (Til), and a 4-day timescale for 
the Fermi-LAT data. The spectral shapes also appear 
different. The AGILE data are characterized by a hard 
curved spectral shape with peak photon energy E p of 
the differential power spectrum E p ~ 300 MeV (Til, 
see Fig. 1). On the contrary, the Fermi-LAT spectrum 
shows a quasi power-law shape extending up to a few 
GeV (All). Without additional analysis, it is not clear 
whether the two datasets are consistent with each other. 
In any case, the hardness of the gamma-ray emission 
and the rapid spectral evolut ion challenge the idealized 
scenario underlying Eq. [1] (Ide Jager fc Hardin g 19921 : 
Ide Jager et al. 19961 lAtovan fc Aharonian 19961) . 

The goal of our paper is twofold: (1) investigate the 
consistency of the published AGILE and Fermi-LAT 
spectral data of the Sept. 2010 event by integrating the 
AGILE data over a 4-day timescale; (2) study the spec- 
tral evolution of a class of synchrotron emission models 
based on freshly accelerated particles in the inner Neb- 
ula, and check its validity for both the AGILE and Fermi- 
LAT data. 

2. SPECTRAL DATA ANALYSIS 

In order to test whether the AGILE and Fermi-LAT 
spectra of the Sept. 2010 event are consistent with each 



other it is necessary to consider data with the same in- 
tegration timescales. In the absence Fermi-LAT spectral 
data on a 2-day timescale, we revisited our AGILE data 
and carried out a 4-day integration which overlaps with 
the Fermi-LAT interval. 

The AGILE 4-day Crab spectrum in the energy range 
50 MeV-3 GeV was obtained by integrating between 
MJD 55457.38-55461.55. We obtained the nebular 
contribution by subtracting from the total emission the 
pulsar contribution corresponding to a flux F(E > 100 
MeV)= 210 ± 30phcm _2 s _1 characterized by a power- 
law spectrum with photon index T = 2.0 ± 0.1 in the 
energy range 50 MeV - 3 GeV. 

Fig. [5] shows the result of our additional analysis of 
the AGILE data together with the published Fermi-LAT 
data. The two data sets are now temporally homoge- 
neous and appear to be in agreement within the errors. 
The apparent power-law behavior of the 4-day spectrum 
in the energy range 50 MeV - 2 GeV can be explained by 
a fast- rise-synchrotron-cooling model (see below). The 
solid curve of Fig. [5] shows the result of our modelling 
for a 4-day average of the rapidly varying spectrum. 

3. A FAST-RISE SYNCHROTRON COOLING MODEL 

We assume that a fresh population of impulsively accel- 
erated electrons/positrons is produced in the inner neb- 
ula within a timescale short compared with all other rel- 
evant cooling timescales. The model presented here has 
general validity and does not depend on a specific site 
in the Nebula as long as the general characteristics of 
the emission fit our assumptions. We assume an efficient 
particle acceleration mechanism 13 that applies simulta- 
neously in one or more contiguous nebular sites that are 
subject to plasma instabilities and/or substantial pul- 
sar wind particle density enhancements. A fraction of 
the total volume of the inner nebula is affected by the 
flaring instability. Consequently, only a fraction of the 
total number of radiating nebular particles contributes 
to the flare. This fact follows from the observed short 
flaring rise-time. For simplicity, we assume a Doppler 
factor V = (1 — f3 cos d)^ 1 ~ 1; a larger Doppler fac- 
tor would imply a smaller particle number N e oc 2?~ 3 , 
a larger emitting region L cx X>, and smaller rest-frame 
particle energies 7 oc T)^ 1 / 2 . 

In our analysis we considered different values of the lo- 
cal magnetic field Bi oc . The rapid observed cooling (2-3 
days) for reasonable values of particle energies imply that 
the local magnetic field is substantially enhanced com- 
pared with B (Til, All). Reconciling the synchrotron 
cooling timescale r s ~ (8 • 10 8 sec) Bf Q 2 7 _1 (where the 
local magnetic field Bi oc is in Gauss, and 7 is a typical 
particle Lorentz factor) with the Sept. 2010 observations 
implies, for 7 rs 5 x 10 9 of electrons irradiating in the GeV 
range, a local magnetic field Bi oc ~ 10~ 3 G that is ~ 5 
larger than the nebular average. 

Our best modelling assumes an emitting region of size 
L = 7 x 10 15 cm, and an enhanced local magnetic field 
Bioc — 10 -3 G that we keep constant in our calculations. 
The acceleration process produces, within a timescale 
shorter than any other relevant timescale, a particle en- 

13 We leave for other investigations the crucial issue of explaining 
the type of plasma wave turbulence leading to the short accelera- 
tion timescale (1 day or shorter). 
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ergy distribution that we model as a double power-law 
distribution (Til) 

dj (7/7i)« + (7/7fc) P2 

where n is the particle number density. The assumption 
T> ~ 1 together with the constraint on _B/ oc implies a 
break energy 7b around 2 x 10 9 , and a normalization fac- 
tor K around 5 x 10~ 10 cm -3 . If the gamma-ray flare is 
related with the persistent local enhancement detected in 
the anvil region by HST and chandra (Til), we can con- 
strain pi = 2.1 and P2 — 2.7, with the particle Lorentz 
factor 7 ranging from j m i n — 10 6 to 7 mQX = 7 x 10 9 . The 
double power-law distribution of Eq. [5] implies maximal 
synchrotron emission between 7b and J ma x an d the total 
particle number required to explain the flaring episode 
turns out to be iV e _/ e+ = J dV (dn/d-y) d"f = 2 x 10 42 , 
where V is an assumed spherical volume of radius L. 

Based on standard synchrotron emissivity and parti- 
cle cooling, we calculated both the particle distribution 
and the photon spectrum evolution keeping Bi oc = lmG 
constant. We show in Fig. [3] the calculated photon spec- 
tra at four different times corresponding to days 1-2-3- 
4. Given our model parameters, fast spectral evolution 
takes place, and the flaring phenomenon fades away 
within the fourth day. We also calculated time spec- 
tral averages of the differential gamma-ray energy flux 
dF / dE for different integration time durations T accord- 
ing to the formula 

with the particle energy loss rate 7 = — "f/r s where r s is 
the synchrotron cooling time. We use the time-integrated 
spectral function of Eq. |3]to model the 2-day (Fig. [JJ and 
4-day (Fig. [2]) integrated spectral data of AGILE and 
Fermi-LAT. We find that the synchrotron peak photon 
energy during day no. 1 is 

E peak = § H — l 2 max * 800 MeV (4) 
I m e c 

which is in good agreement with the peak shown in the 2- 
day averaged AGILE spectrum (Fig. [JJ). In the absence 
of very strong Doppler effects, our measured spectrum 
and the calculated E pea k violate the expectations from 
Eq. [TJ Doppler effects with V ~ a few would not alter 
this conclusion. We find that the emission from inverse 
Compton scattering of the flaring particle population is 
negligible. 

We also note that the spectral shape calculated in Fig. 
[3]and measured in Fig. Q]contradicts a simple translation 
by a Doppler factor of the average nebular data showing 
the synchrotron burn-off. The additional population of 
energized particles necessary to explain the flare can be 
successfully modelled by Eq. [3] (that can also account 
for the X-ray and optical "afterglow" in the anvil region 
measured by Chandra and HST, Til). Our model ap- 
plies to emitting regions idealized as standing sites or as 
regions within MHD outflows (such as the anvil features). 
Adiabatic expansion could play a role in contributing to 
the gamma-ray flux decrease for an emitting site speed 
of order of the sound speed. We checked the relevance 
of the adiabatic expansion in our model and concluded 



that it could explain a good fraction of the observed flux 
decrease. In this case, our estimate of the the local mag- 
netic field Bi oc would be an upper limit. Moreover, there 
would be no direct relation between the gamma-ray emis- 
sion and the X-ray emission of the anvil enhancement be- 
cause, at variance of the synchrotron cooling, adiabatic 
models cannot account for the persistent brightening in 
the X-rays (see Fig. [3]). 

We note that also a purely Maxwellian particle energy 
distribution (presented in Til, and resulting from par- 
ticle shocks with no non-thermal tails) can also in prin- 
ciple explain the spectral evolution above 100 MeV: also 
in this case there would be no direct spectral connection 
between the gamma-ray emission and the X-ray/optical 
properties of localized regions in the Nebula. 

4. AN EVEN FASTER EVOLUTION MODEL 

In our study we considered also the possibility of a 
flux and spectral evolution even faster than that shown 
in Figs. 1-3. The analysis of the Fermi-LAT Sept. 2010 
data by Balbo et al. 2010 suggests indeed that the spec- 
tral evolution may occur on an overall timescale even 
faster than 1-2 days. A cooling timescale of ~ lday is 
explained in our synchrotron model for Bi oc ~ 2.5 mG, 
and j m ax = 5 10 9 . In this case, the overall flaring episode 
lasting ~ 4 days is characterized by a sequence of short 
acceleration and cooling episodes lasting 1-2 days. Our 
analysis remains valid also in presence of a faster evo- 
lution. In this case, the_ magnetic field is determined to 
have a value Bi oc ~ 10 B. 

We note that slower events lasting 6-7 days, as the 
2007 flare presented in (Til), can be interpreted by the 
same model with larger region involved L w 5 10 16 cm 
and similar Bi oc . Determining the gamma-ray temporal 
structure on timescales shorter than 2 days is limited 
by photon statistics. The time-resolved analysis of the 
AGILE gamma-ray data for the Sept. 2010 event will be 
presented elsewhere. 

5. DISCUSSION AND CONCLUSIONS 

The Sept. 2010 event of the Crab Nebula lasting 
~ 4 days is currently the shortest detected gamma-ray 
flare. Our analysis shows that the flux and spectral evo- 
lution of this event are well described by a model charac- 
terized by very fast (shorter than ~ lday) particle accel- 
eration and by synchrotron cooling in a local magnetic 
field 5-10 times larger than the average nebular value 
B. Both the AGILE and Fermi-LAT gamma-ray spec- 
tral data are consistent with each other within a 4-day 
timescale. Our analysis of the AGILE data on a 2-day 
timescale clearly shows that the emission is peaked at 
the photon energy of Eq. 21 which is almost one order of 
magnitude larger than the "synchrotron burn-off" con- 
straint of Eq. [TJ The flaring mechanism in the Crab 
Nebula is quite remarkable: it accelerates particles to 
the largest kinetic energies (PeV) associable to a specific 
astrophysical source and does it within the shortest time 
ever detected in a nebular environment. 

Our results challenge the physical assumptions un- 
derlying Eq. [JJ and in particular acceleration models 
based on "slow" processes. As we showed above, ex- 
planations in terms of Doppler boosting are problem- 
atic in light of the measured spectral curvature of the 
AGILE data. Even though a theoretical study of pos- 
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sible acceleration mechanisms consistent with the data 
discussed here is beyond the scope of this paper, we 
can briefly mention some of the difficulties. First-order 
Fermi acceleration with particles gaining energy by dif- 
fusing stochastically back and forth a shock front (e.g., 
Blandford & Ostriker 1978, Bell 1978, Drury 1983) ap- 
pears to be too slow and is drastically challenged by 
our findings. In particular, it is difficult to see how a 
diffusive shock acceleration mechanism can violate Eq. 
[TJ A locally enhanced (over Bi oc ) electric field can pro- 
duce a sort of "runaway" of kinetic energy gains with 
an acceleration rate larger than the synchrotron cool- 
ing rate. However, despite some attempts and analogies 
with other astrophysical contexts (e.g., pulsar magneto- 
spheres), it is currently not clear how this mechanism 
can be implemented in the Crab Nebula. MHD mod- 
els of the pulsar wind (e.g., Komissarov & Lyubarsky 
2004, Del Zanna et al. 2004, Camus et al. 2009, Komis- 
sarov & Lyutikov 2010), address the turbulence and the 
limit-cycle behavior of the instabilities. These features 
may in principle favor substantial local magnetic field en- 
hancements. However, the calculated timescales of these 
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instabilities (e.g., Camus et al. 2009) are several or- 
ders of magnitudes longer than what we detected in the 
Crab Nebula. Shock-drift acceleration (|Kirk et al. 2000) 
tends to occur on a timescale shorter than for diffusive 
processes. However, it is not clear whether the required 
efficiency can be reached in the flaring Crab Nebula site, 
and whether Eq. @] can be obtained. Shocks medi- 
ated by ions in the pulsar wind that resonantly acceler- 
ate pairs by magnetosonic w aves (iGallant fe Arons 19"95 
ISpitkovskv fe Arons 2004 lArons 20081 ) are typically 
slow, and are most likely not applicable in the X-ray 
and optically enhanced pulsar polar jet regions of Til. 

The challenge provided by the Crab Nebula gamma-ray 
flaring requires a thorough investigation of the mecha- 
nisms leading to efficient particle acceleration and to a 
natural justification of Eq. |4] The issue will be eluci- 
dated by future Chandra X-ray and HST optical obser- 
vations of the inner Crab Nebula that will be carried out 
in search of the gamma-ray flaring site. 
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Fig. 1. — AGILE (filled squares) 2-day averaged data of the Sept. 2010 gamma-ray flare of the Crab Nebula. Pulsar data have been 
subtracted. The solid line represents the 2-day averaged synchrotron emission model of the Sept. 2010 flare summed with the standard 
nebular emission as discussed in the text. The dashed curve in red shows the flaring component averaged over 2 days. Data points in open 
circles give the standard average Crab Nebula spectrum that we model by the dashed black curve. The spectral region marked in green 
shows the X-ray data of "source A" of Til. 
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Fig. 2. — AGILE (filled squares) and Fermi-LAT (open triangles) 4-day averaged spectral data of the Sept. 2010 gamma-ray flare of 
the Crab Nebula. Pulsar data have been subtracted. The solid line represents the 4-day averaged synchrotron emission model of the Sept. 
2010 flare summed with the standard nebular emission as discussed in the text. The dashed curve in green shows the flaring component 
averaged over 4 days. Data points in open circles give the standard average Crab Nebula spectrum that we model by the dashed black 
curve. 
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Fig. 3. — Spectral evolution of the gamma-ray Crab Nebula September 2010 flare as obtained by our fast-rise-synchrotron-cooling model. 
The upper curve shows the spectrum at the starting time, and the lower curves show the spectra after 1, 2, 3, and 60 days respectively. 
Note the persistence of the X-ray emission in the "source A" of Til localized by Chandra. 



